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We develop a new approach to the theoretical treatment of the separatrix 
chaos, using a special analysis of the separatrix map. The approach allows us 
to describe boundaries of the separatrix chaotic layer in the Poincare section 
and transport within the layer. We show that the maximum which the width 
of the layer in energy takes as the perturbation frequency varies is much larger 
than the perturbation amplitude, in contrast to predictions by earlier theories 
suggesting that the maximum width is of the order of the amplitude. The 
approach has also allowed us to develop the self-consistent theory of the earlier 
discovered (PRL 90, 174101 (2003)) drastic facilitation of the onset of global 
chaos between adjacent separatrices. Simulations agree with the theory. 
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1. Introduction 

Even a weak perturbation of an intcgrable system possessing a scparatrix 
results in the onset of chaotic motion inside a layer which we shall fur- 
ther call as the separatrix chaotic layer (SL). The separatrix chaos plays 
a fundamental role for the Hamiltonian chaos, being also relevant to vari- 
ous applications.^"^ The boundaries of the SL in the Poincare section may 
be easily found numerically."^ However, it is also important, both from the 
theoretical and practical points of view, to be able to theoretically calculate 
them and describe transport within the SL. 

One of the most powerful theoretical tools for the SL study is the sepa- 
ratrix map (SM), introduced in^ for the nearly intcgrable systems with the 
3/2 degrees of freedom and called sometimes^ as the Zaslavsky separatrix 
map. It may also be generalized for systems with more degrees of freedom 
and for strongly non-integrable systems (see^ for the most recent major 
review). We shall further consider the case of the 3/2 degrees of freedom 
but the generalization of our method for other cases may be done too. 

One of the most interesting for physical applications relevant quantities 
is the width of the SL in energy. ^""'■^"^^ There are various heuristic crite- 
ria^"'* based on the separatrix map and various conjectures. The width by 
these criteria does not depend on the angle and, as a function of a per- 
turbation frequency cof, possesses a maximum at cof of the order of the 
eigenfrequency in the stable state uq while the maximum itself is of the or- 
der of the perturbation amplitude h. However, the work*° has demonstrated 
in simulations for double-separatrix systems that the maximum width may 
be much larger as the SL absorbs one or two nonlinear resonances. The 
recent work^^ has proved this, developing a new method for the analysis of 
the separatrix map. The method is of a general validity, as shown in the 
present work. We show that the maximum width occurs at the frequency 
which is typically smaller than uiq by the logarithmic factor ln(l//i) while 
the maximum width is typically much larger than h - either by a numerical 
factor or by the logarithmic factor (apart from the adiabatic divergence in 
certain class of systems^). Besides, the method allows to describe major 
statistical properties of transport within the SL. 

Note that there were various mathematical works considering the SL in 
rather different contexts (see^ for the review). In particular, they analysed 
the SL width in normal coordinates. However, to the best of our knowl- 
edge, these works do not specify the relation between the normal coordi- 
nates and variables conventional in physics (e.g. energy-angle or coordinate- 
momentum). Besides, these works just estimate the width from above and 
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below while our method allows to carry out an accurate calciilation of the 
width in energy and, moreover, of the SL boundaries in the Poincare sec- 
tion. Finally and most importantly, the methods described in^ do not make 
a resolution between the resonance frequency range and other frequency 
ranges while our method shows that the SLs in these ranges drastically 
differ from each other, as confirmed by simulations. 

Below, we describe the basic ideas of our method (Sec. 2), review the 
results of its application to the doublc-scparatrix case (Sec. 3) and present 
rough estimates for the single-separatrix case (Sec. 4). 



2. Basic ideas 

Consider any ID Hamiltonian system possessing at least one separatrix. 
Let us add a weak time-periodic perturbation. 



H = Ho{p,q) + hV{p,q,t), V {j>, q,t + 2tt / f) = V{p,q,t), /i< 1. (1) 

The motion near any of the separatrices may be approximated by the 
separatrix map (SM).!"*^'!^ The map slightly differs for different types of 
separatrix. Our method applies to all types but, to be concrete, we consider 
in this section only the separatrix with a single saddle and two loops (like 
in a double- well potential system). Then the SM reads as^^ (cf. also^"^): 



Ei+i = Ei + aih€sm{(pi), (2) 
w/7r(3 - sign(£;;+i - Eg)) 

= + M^) ' 

CTj+i = cr, sign(£;s - Ei+i), 

/OO 
dt dHo/dp\j^^sm{ujft), 
-co 

Ei = Ha(p,q)\^^ , = u>fti, Ui = sign(ai/o/5p|tJ, 

where Es is the separatrix energy. The quantity e is often called as the Mel- 
nikov^ or Poincare-Melnikov^ integral. The quantity ^ = /i|e| is sometimes 
called separatrix split.^ For the sake of simplicity, let absolute values of all 
parameters of Hq and of be '--^ 1. Then |e| ^ 1 too, if w/ ~ 1. 
Consider two most general ideas of the method. 

1. The SM trajectory that includes any state with E = Eg is chaotic 

since the angle of this state is not correlated with the angle of the state at 
the preceding step of the map, due to the divergence of lv~^{E Eg). 
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2. The frequency of eigenoscillation as a function of energy is propor- 
tional to the reciprocal of the logarithmic factor: 



^(^) = mW^^' " = ^ ' 

\E-E,\ « Aif ~ - ~ - 

fx 2) 

{Elf ' are energies of the stable states) . 

Given that the argument of the logarithm is large in the relevant range 
of E, the function lo{E) is nearly constant at a rather significant variation 
of the argument. Therefore, as the SM maps the state {Eq = Eg, <po) onto 
the state with E = Ei = Eg + (To/ie sin((/)o), the value of uj{E) for the given 
sign(c7oesin((/?o)) is nearly the same for most of angles ipo (except the close 
vicinity of tt multiples), namely 



uj{E) K, a;(='=) = uj{Es ± h) for sign((Toesin((^o)) = ±1- (4) 

Moreover, even if the deviation of the trajectory of the SM from the sep- 
aratrix further increases /decreases, uj{E) remains close to Wr^^ provided the 
deviation is not too large/small, namely if | h\{{E — E g) / h)\ <C ln(A_H'/ft,). If 
w/ ~ uJr'^\ then the evolution of the map (2) may be regular-like for a long 
time until the energy returns back to the close vicinity of the separatrix, 
where the correlation of angle is lost again. Such a behavior is especially 
pronounced if the perturbation frequency is close to u}\^'^ or uji ^ or to one 
of their multiples of a not too high order: the resonance between the per- 
turbation and eigenoscillation gives rise to an accumulation of the energy 
gain for many steps of the SM, which results in the deviation of E from Eg 
greatly exceeding the separatrix split (5 ~ /i. As a function of w/, the largest 
(along the SL boundary) deviation from the separatrix takes its maximum 
at the frequencies slightly exceeding Wr^^ and ' , for the upper and lower 
boundaries of the SL respectively. This corresponds to the absorption of the 
nonlinear resonance by the SL. 

The description of the regular-like parts of the chaotic trajectory in the 
case close to the resonance may be done within the resonance approxima- 
tion. The explicit matching between the SM and the resonance approxi- 
mation is carried out in.^^ The resonance approximation is done in terms 
of slow variables, action I = I{E) (note that dl/dE = co~^{E)) and slow 
angle ip = — Wft = ip — (p, by means of the resonance Hamiltonian 
F(J,^/').i-4.9-ii 
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angle \j/ 

Fig. 1. Schematic illustrations to the formation of the SL upper boundary. The sepa- 
ratrix energy Es corresponds to the lower boundary of the box. The GSS curve is shown 
by the dashed line. Solid lines show examples of resonant trajectories (RTs) overlapping 
the GSS curve. The SL boundary (thick solid line) is formed by: (a) the RT tangent to 
the GSS curve, or (b) the upper part of the self-intersecting RT (resonant separatrix). 



Fig. 1 schematically illustrates the formation of the upper boundary of 
the SL in the Poincarc section presented in the E — ip variables. The chaotic 
trajectory jumps from the separatrix onto the curve which we call the upper 
generalized separatrix split (GSS) curve, within an even tt interval: 



Eglii') = E, + d\sin{^)\, S = h\e\, (5) 
Vi e [7r + 27rn,27r + 27rn], n = 0, ±1, ±2, ... 

The GSS curve relates to the SM equation for E with Eq = Eg. The rel- 
evance of just the even tt intervals of is a consequence of a necessarily 
positive sign of Ei — Eq as far as we consider the upper boundary of the 
SL. Then the system follows the trajectory of the resonance Hamiltonian H 
{E necessarily increases initially because V'(O) is within an even tt interval), 
until it again reaches the GSS curve (necessarily in an odd tt interval, where 
E decreases). After that, the system jumps onto the separatrix, where the 
angle correlation is lost. There is a statistical distribution of regular-like 
parts of the chaotic trajectory corresponding to the homogeneous distribu- 
tion of initial angles. The boundary of the SL is formed by the trajectory of 
the resonance Hamiltonian (RT) starting from the GSS curve at such ini- 
tial angle that provides for the deviation of the RT from the separatrix to 
be larger than that for any other initial angle. From the topological point 
of view, there may be two different situations: either the relevant RT is 
tangent to the GSS curve (Fig. 1(a)) or it intersects it while being a self- 
intersecting trajectory i.e. a separatrix of the resonance Hamiltonian (Fig. 
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1(b)). In the latter case, the outer boundary is formed by the outer part of 
the separatrix (in Fig. 1(b), relating to the upper boundary, this part lies 
above the saddle s; an example of the boundary formed by the separatrix 
of a different type is shown in Fig. 4(a)). 

The boundary strongly depends on the angle. The maximal/minimal 
deviation of E from Eg is much larger/smaller than S. 
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Fig. 2. The potential U (q) = (0.2 — sin(g))^/2, the separatrices and the eigenfrequency 
U!{E) of the unperturbed system Ho = p^/2 -f C/(q), in (a), (b) and (c), respectively. 
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Fig. 3. (a). The diagram indicating (shading) the perturbation parameters range for 
which global chaos exists in the perturbed system, H = Ho(p,q) — hq cos{lu ft), (b). The 
comparison of the major minimum of the diagram with the lowest-order theory (dashed 
lines) and the theory allowing for higher-order corrections (solid lines). 



3. Application to the double-separatrix chaos 

It has been found in^'' that, if the unperturbed Hamiltonian Hq possesses 
more than one separatrix (cf. Fig. 2) while the perturbation is time-periodic, 
the onset of global chaos in between the adjacent separatrices possesses a 
remarkable feature: it is greatly facilitated if the perturbation frequency 
w/ is close to certain frequencies: the perturbation amplitude h required 
for the chaos onset is much smaller for such frequencies than for others 
(Fig. 3). This is related to the characteristic shape of oj{E) in between the 
separatrices: Cij{E) approaches a rectangular ioira in the asymptotic limit of 
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a small separation between the separatrices (cf. Fig. 2(c)). If cuf is slightly 
smaller than the local maximum Um of lo{E), then there are two nonlinear 
resonances that are very wide in energy: they may simultaneously overlap 
each other and the separatrix chaotic layers, that occurs at the value of h 
which is logarithmically smaller than a typical value required for the chaos 
onset when ujf lies beyond the immediate vicinity of ujm- 

Using the semi-heuristic approach, the lowest-order asymptotic theory 
based on the resonance Hamiltonian analysis was developed in^*^ for the 
minima of the global chaos boundary hgdujf). The theory was compared 
to results of computer simulations for the given example (for a moderately 
small separation: Fig. 2). The value of in the minimum was well described 
by the lowest-order formula but the descrepancy for hgc in the minimum 
was nearly 50%. Besides, it was unclear how the overlap of resonances with 
the chaotic layers occured and why even a small excess of h over hgc resulted 
in the onset of chaos in a large area of the Poincare section despite that 
chaotic layers associated with the nonlinear resonances were exponentially 
narrow for h < hgc- These problems have been resolved in our recent work^^ 
developing the method similar to that described in the previous section. The 
agreement between the theory and simulations has greatly improved (Fig. 
3(b)). 

Our present work generalizes the above method for any separatrix (see 
Sec. 2 above and Sec. 4 below). The general validity of the method is 
brightly demonstrated by Fig. 4, that shows the direct comparison of the 
theoretical SL boundaries and the SLs generated by computer. Though the 
layers are still related to the system with two separatrices, the given per- 
turbation amplitude is so small that the layers are well separated from each 
other and therefore the presence of a second separatrix does not play a sig- 
nificant role for any of the layers. The lower and upper layers demonstrate 
two characteristic situations: the lower layer has not yet absorbed the rele- 
vant nonlinear resonance (though the closeness to the resonance gives rise 
to a rather strong increase of the layer width) while the upper layer has 
absorbed the resonance so that its maximal width greatly exceeds that of 
the lower layer and there is a large island of stability in the layer. 

4. Single-separatrix layer: estimates of the largest width 

As already mentioned, the SL width in energy takes its maximum at w/ w 
uj'^^~\ The rigorous treatment for the single-separatrix case may be done 
similar to the double-separatrix case.^'^ It will be done elsewhere while, 
here, we present rough estimates for the width. 
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Fig. 4. The separatrix chaotic layers (shaded) in the plane of a<;tion I and slow angle ip 
for the system exploited in Fig. 3, for h = 0.003 and ujf = 0.401, as described by our the- 
ory. The dashed lines represent the relevant GSS curves. The solid lines represent relevant 
trajectories of the resonance Hamiltonian H: two solid lines with the saddles represent 
separatrices of the nonlinear resonances, two other solid lines are the trajectories of H 
which are tangent to the lower and upper GSS curves respectively, (b). Comparison of 
the layers obtained from computer simulations (dots) with the theoretically calculated 
boundaries (solid lines) shown in the box (a). 

Let US transform from p — q to action / = I{E) and angle and expand 
V in Eq. (1) into the double Fourier series (in t and tp)-} 



Let us single out the term with the maximum absolute value of Vk,i 
(typically, it corresponds to k = I = 1) and denote it as Voe'^: 



The maximum width of the SL corresponds to the perturbation fre- 
quency at which the SL has just absorbed the widest nonlinear resonance 
(cf.^^). From the rigorous results for the double-separatrix case,^-'^ we may 
assume that the SL width is dominated by the width of the nonlinear reso- 
nance then, i.e. the width of the resonance separatrix in energy AE greatly 
exceeds the minimal separation in energy between the resonance separatrix 
and Eg. Obviously, this assumption should be verified in the end. 

Strictly speaking, Vb strongly depends on / in the relevant range of 
/, end the rigorous treatment of the nonlinear resonance is complicated 
(cf.^^^^). But for a rough estimate of the resonance width, it is suffi- 
cient to use a simple Chirikov approximation of the resonance Hamiltonian 





k,l 



Vo{I) = max{\Vk,i\). 

k,l 



(7) 
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H{I,tjj),^ 4,9,12 -^iiicii reduces to the auxiliary pendulum dynamics. The 
width of the corresponding resonance separatrix in energy is expressed as 



The Chirkov approximation assumes that |dw/d£^| « const within the 
range of energies relevant to the resonance separatrix. In our case, it is not 
so since the quantity jdw/di^l w uj'j- / {aTruJo\E — Eg]) strongly varies within 
the relevant range of energies. However, we may still use Eq. (8) for the 
rough estimate, putting in it \E — Eg] ~ AE. Then we obtain for AE the 
following rough asymptotic equation 



AErs.Vo{E = E,±AE)hln{l/h), h^O, (9) 

where we took into account that the relevant w/ is close to 
LOo/ln{l/h) and omitted numerical factors. 

The asymptotic solution of Eq. (9) essentially depends on the function 
Vo{x). In this context, all perturbed systems may be divided in two classes. 

1. The separatrix of the unperturbed system has more than one saddle 
while the relevant coefficient Vi = Vi {E, ^) in the Fourier expansion (in 
time) of the perturbation V possesses different values in adjacent saddles. 
An archetypal example is a pendulum subject to a dipole time-periodic per- 
turbation.^'^ If E is close to Eg, then the system stays mostly near one of 
the saddles, so that Vi depends on ijj nearly in a piece- wise manner: it oscil- 
lates between the values corresponding to the adjacent saddles. Therefore, 
Vq (which is the absolute value of the relevant coefficient in the Fourier ex- 
pansion of Vi in ip) approaches in the asymptotic limit h ^ some non-zero 
constant. As follows from Eq. (8), AE' is logarithmically large: 



AE ^h\n{l/h)->h, h ^ 0. (10) 

This estimate agrees with the rigorous result and the result of simulations 

in the case considered in.^-'^ 

2. Either the separatrix has a single saddle (like for a double- well po- 
tential system^'*) or the separatrix has more than one saddle while the per- 
turbation possesses identically equal values at different saddles. Archetypal 
examples arc a pendulum in a wave with the same spatial period"'^^* and 
a pendulum with the oscillating suspension point. ^ Then Vi{E — > Es,tp) 
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keeps nearly one and the same value for most of the period of ip (as it 
stays most of the period near the saddle/s): it significantly differs from this 
value only during a small part of the period, which is uj{E)/loq. Hence, 
Vo{E = Es± AE) ~ 1/ ln(l/Ai;), so that the solution of Eq. (9) is: 

AE ^ h, (11) 

This means that the asymptotic functional dependence of the resonance 
width is the same as that of the SL width in frequency ranges beyond 
the resonance. So, the functional dependence of the SL in the resonance 
range is of the same type as beyond it, being oc h. At the same time, a 
ratio between the resonance width and h may still be a large number. Both 
these conclusions are in agreement with computer simulations. Thus, for 
the archetypal example of the Duffing oscillator subject to the dipolc time- 
periodic perturbation,^'* the ratio AE/h approaches in the limit h ^ the 
constant value approximately equal to 20: see Fig. 3(b) in.® For another 
archetypal example, namely a pendulum in the wave of the same spatial 
period, our recent simulations for the parameters exploited in^'^''^ have 
shown that AE/h '^—^ const « 50. 
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